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A year after Fermi was launched, the number of known gamma-ray pulsars has increased dramatically. For 
the first time, a sizable population of pulsars has been discovered in gamma-ray data alone. For the first time, 
millisecond pulsars have been confirmed as powerful sources of gamma-ray emission, and a whole population 
of these objects is seen with the LAT. The remaining gamma-ray pulsars are young pulsars, discovered via 
an efficient collaboration with radio and X-ray telescopes. It is now clear that a large fraction of the nearby 
energetic pulsars are gamma-ray emitters, whose luminosity grows with the spin-down energy loss rate. Many 
previously unidentified EGRET sources turn out to be pulsars. Many of the detected pulsars are found to 
be powering pulsar wind nebulae, and some are associated with TeV sources. The Fermi LAT is expected 
to detect more pulsars in gamma rays in the coming years, while multi-wavelength follow ups should detect 
Fermj-discovered pulsars. The data already revealed that gamma-ray pulsars generally emit fan-like beams 
sweeping over a large fraction of the sky and produced in the outer magnetosphere. 



1. Introduction 

After the Compton Gamma Ray Observatory 
(CGRO) ended its activity in 2000, seven rotation- 
powered pulsars were known to emit in gamma rays, 
all normal pulsars, plus a couple of low-significance 
detections, inclu ding the millisecond pulsar PSR 
J0218+4232 (see ^ for a review of CGRO results 
on pulsars). More recently, the AGILE telescope [11] 
reported the detection of two other normal pulsars in 
gamma rays [13, HI] plus a few other marginal de- 
tections, including that of another millisecond pulsar 
[33| . yet to be confirmed. Less than ten pulsars were 
hence known to emit in gamma rays at the beginning 
of the Fermi mission, a small number compared to the 
almost 2000 known at radio wavelengths. 

Because of this limited sample, all attempts to 
model the fundamental phenomena producing high- 
energy emission in pulsar magnetospheres were so far 
confronted by a severe lack of observational data. 
Namely, it was unclear whether high-energy photons 
are produced above the magnetic poles, or rather in 
gaps formed at higher altitude in the magnetosphere 
(see (28j for a review of pulsar gamma-ray emission 
models). These models predict different features in 
pulsar light curves (signal intensity as a function of 
the rotational phase): number of gamma-ray peaks 
if any, phase separation between peaks if more than 
one, and phase separation between components at dif- 
ferent wavelengths. Emission spectra also differ from 
one model to one another: for normal pulsars emission 
close to the surface would suffer magnetic pair atten- 
uation, leading to a super-exponential energy cutoff 
in spectra, while for higher altitude emission the ex- 
pected energy cutoff would be simple-exponential. Be- 
cause they assume different emission geometries, the 
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models predict different ratios of radio-loud and radio- 
silent gamma-ray pulsars, as their radio beam may 
or may not sweep across the Earth. Before Fermi, 
only one such radio-silent gamma-ray pulsar had been 
detected, the Geminga pulsar, so that this ratio re- 
mained largely unconstrained. Finally, it was unclear 
whether millisecond pulsars (MSPs) could emit high- 
energy emission, and if they do, how their emission 
would compare to that of normal gamma-ray pul- 
sars. The prospects for detections of MSPs were high 
though, after the few marginal detections in EGRET 
and AGILE data. The launch of the Fermi obser- 
vatory on 11 June 2008 provided a new opportunity 
to study the high-energy emission from pulsars and 
address the above questions. 

The Large Area Telescope (LAT) :20j , main instru- 
ment of the Fermi observatory, is sensitive to pho- 
tons between 20 MeV and more than 300 GeV. It has 
an unprecedented sensitivity, thanks to its large field 
of view of ^ 2.4 sr, large effective area (~ 8000 cm^ 
above 1 GeV on axis) and improved angular resolution 
(0.6° at 1 GeV). It also has an excellent timing accu- 
racy, better than 1 iis [ll| , crucial for pulsar observa- 
tions. Furthermore, the survey observation strategy 
of the Fermi observatory makes the LAT a prime in- 
strument for the discovery of new sources. We review 
here the results of pulsar observations with the Fermi 
LAT during its first year of activity, and discuss the 
vast prospects following these results. 



2. Fermi LAT observations of pulsars 
2.1. Observations of known pulsars 

There are two different approaches for the obser- 
vation of pulsars with Fermi. One is to search for 
gamma-ray pulsations from known pulsars, by phase- 
folding the photon arrival dates using pulsar rotation 
parameters measured at other wavelengths, namely 
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Figure 1: Spectral energy distribution for the PSR 
J2021+3651 pulsar. The solid line shows the total flux 
emitted by the pulsar. The dotted and dot-dash lines 
correspond to the first and second peak in the gamma-ray 
light curve. See d for additional details on this Figure. 



radio and X-rays. Because radio and X-ray telescopes 
cannot monitor all ~ 2000 known pulsars, a list of ^ 
230 prime candidates for detection in gamma rays has 
been defined prior to the launch, including the seven 
CGRO detections [32] . These prime targets were then 
monitored by large radio telescopes around the world, 
as well as X-ray telescopes, aiming at measuring their 
rotation parameters (pulsar ephemerides) and keeping 
them accurate over the Fermi observation range. Pa- 
rameters were also obtained for a priori secondary tar- 
gets, expanding the number of pulsars to be searched 
for pulsations with the LAT to more than 700 objects. 

Using these parameters, the LAT detected gamma- 
ray pulsations from a number of normal radio-loud 
pulsars, all highly energetic. Two of these pulsars, 
J0659-hl414 and J1048-5832 [10], had been seen 
with EGRET with low significance. Other objects lie 
in the error box of EGRET unidentified sources [2^ : 
J1028-5819 in 3EG J1027-5817 [|, J1420-6048 in 
3EG J1420-6038 13, J2021-h3651, also seen with the 
AGILE telescope Till, in 3EG J2021-h3716 [| and 
J2229-h6114, in 3EG J2227-I-6122 [l^. The Fermi 
LAT has so far detected a total of 22 of these ob- 
jects with high-confidence, over 100 MeV Table 
3 of [13] gives the gamma-ray pulse shape parameters 
for these pulsars. Although the nature of a few light 
curves are unclear because of limited photon statistics, 
the typical light curve is two-peaked, with a separation 
between these peaks of ^ 0.4, and the first gamma- 
ray peak lagging the main radio component by 0.1 to 
0.2. Good examples of this typical shape are e.g. the 
Vela pulsar Q or PSR J2021-I-3651 0. A few excep- 
tions to this trend exist though, like the one-peaked 
PSR J2229-H6114 [13], for which the unique gamma- 
ray peak lags the radio component by approximately 
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Figure 2: Exponential cutoff energy as a function of the 
rotational phase, for the Vela pulsar. Figure adapted 
from ,40i]. 



half a rotation. The emission spectra have been ex- 
amined as well (see Table 5 of [IJ] for a summary of 
spectral results). As illustrated in Figure [T] for PSR 
J2021-t-3651, the energy spectra are found to be gen- 
erally well fit by exponentially cutoff power laws, with 
cutoff energies below 10 GeV, in agreement with what 
had been measured with the EGRET telescope [39| . 

Among radio-loud normal pulsars, the ones de- 
tected by the EGRET telescope or other previous 
experiments are in general prime targets for detailed 
spectral analyses because of their brightness, although 
the LAT has seen other bright pulsars that EGRET 
would have detected if accurate rotational parameters 
were available at that time, such as PSR J2021-I-3651. 
Early in the Fermi mission, the high-resolution light 
curves produced with LAT data for the Vela pulsar 
revealed a small third peak between the two main 
gamma-ray peaks, shifting to larger phases with in- 
creasing energy [6| . The large amount of gamma-ray 
photons accumulated during this first year of data tak- 
ing now enables very fine analysis of its emission prop- 
erties as a function of rotational phase. Figure [2] shows 
the energy cutoff of the emitted spectra as a function 
of phase. The important variations of the cutoff en- 
ergy between the two main peaks is consistent with 
the above description of a third peak shifting towards 
higher rotational phase with increasing energy. Simi- 
larly complex behaviors of the spectral index and the 
cutoff energy are observed for other bright EGRET 
pulsars such as Geminga, J1057-5226, J1709-4429, 
or J1952-h3252 [11], providing valuable input for 
the theoretical models. 

The search for pulsed gamma-ray emission from 
MSPs has been fruitful as well: pulsations were found 
for the nearby isolated MSP PSR J0030-f0451 after 
three months of data taking, making this pulsar the 
first MSP firmly detected in gamma rays 0] . Similarly 
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Figure 3: Radio (lower panels) and gamma-ray (upper panels) pulse profiles for three of the nine millisecond pulsars 
detected so far with the Fermi LAT. Two rotations are shown for clarity. Additional details on the construction of 
these light curves can be found in [J. 



to normal gamma-ray pulsars detected with the LAT, 
J0030-K0451 hes in the error box of EGR J0028+0457, 
listed in the revised EGRET catalog of gamma-ray 
sources [2^. After nine months of observation, the 
LAT detected pulsations from eight galactic MSPs, 
including the confirmation of PSR J0218-f4232 as a 
gamma-ray pulsar [1, [3l| . Figure [3] shows radio and 
gamma-ray light curves for a selection of the eight de- 
tected MSPs. The light curves seen so far for MSPs 
resemble those of normal pulsars, as illustrated in Fig- 
ure[3l PSR J0030-I-0451 has two peaks in gamma rays, 
separated by ~ 0.45 and its first gamma-ray peak 
lags the main radio peak by 0.15, making its pul- 
sar shape comparable to that of Vela or J2021+3651. 
PSR J0437— 4715 shows one broad gamma ray peak, 
lagging the radio peak by 0.4, similar to what is 
observed for the normal pulsar J2229-I-6114. PSR 
J2124— 3358 seems to exhibit a more complex light 
curve, as does the radio emission, with a main broad 
peak and possible secondary features. More data will 
help determine the peak multiplicity as is also the 
case for the normal pulsar J1420— 6048 [14|. Sim- 
ilarly, their spectral properties do not differ funda- 
mentally from that of normal pulsars (see Table 1 of 
Q); although super-exponentially cutoffs could not be 
tested, because of the limited size of the gamma-ray 
sample. As the LAT continues to accumulate pho- 
tons, this newly revealed population of gamma-ray 
MSPs grows. For instance, the LAT recently discov- 
ered faint pulsed emission from the binary MSP PSR 
J0034-0534 0. 



2.2. Detections of new pulsars 

The other approach is to search the gamma-ray data 
for yet unknown pulsars. The gamma rays are rare 



though: the LAT typically gets a few hundred pho- 
tons for several months of data taking on a particular 
source. Because of these long integration times, di- 
rect Fourier transforms of the signal are computation- 
ally intensive. One can however search for gamma- 
ray pulsations in differences between photon arrival 
times. This so-called "time-differencing" technique 
dramatically reduces the needed computational power 
[iqI . [isj . Although this technique did not foster new 
detections of gamma-ray pulsars in the EGRET data, 
it has been extremely eflticient with Fermi data: 16 
slowly-rotating pulsars were found after five months 
of gamma-ray observation, including a pulsar in the 
young supernova remnant CTAl, another pulsar coin- 
cident with the isolated neutron star RX J1836-I-5925, 
and other objects likely associated with supernova 
remnants or pulsar wind nebulae [J 0|- In total, 13 
of these 16 LAT-detected pulsars lie in the error box 
of third EGRET catalogue unidentified sources. Fig- 
ured] shows a sample of the observed gamma-ray light 
curves for these 16 pulsars. As expected, these light 
curves resemble those of other gamma-ray pulsars: 
they are generally two-peaked, with a few showing one 
broad peak. Their spectral properties are also consis- 
tent with that of the gamma-ray population (see Ta- 
ble 5 of 14]). More recently, a search for pulsars in 10 
months of data using ~650 LAT source positions unas- 
sociated with possible active galactic nuclei yielded 
the detection of eight new objects [3], expanding the 
number of LAT detected pulsars to 24 objects. Among 
these eight new pulsars, only one might be associated 
with an EGRET unidentified source, indicating that 
the Fermi LAT is perhaps finishing the identification 
of unknown pulsars hidden behind EGRET sources, 
after one year of observation. 

As discussed above, the ratio of radio-loud to radio- 
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Figure 4: Gamma-ray pulse profiles for a subset of the pulsars discovered with the Fermi LAT. Two pulsar rotations 
are shown for clarity. Details on the construction of these profiles can be found in 0]. 
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Figure 5: Gamma-ray light curve of PSR J2032-I-4127 
along with the radio profile measured with the GBT 
radio telescope [S^]. Two other LAT-detected pulsars 
have so far been detected at radio wavelengths, at the 
Green Bank and Arecibo telescopes. 



again very similar to that of other radio-loud pulsars. 
The other GBT-detected pulsar, PSR J1741-2054, is 
faint and close, and may have the lowest radio lu- 
minosity of all known radio pulsars. Similarly, the 
J1907-I-0602 detected with the Arecibo telescope is 
extremely faint 4 /iJy flux density at 1.4 GHz), al- 
though the distance inferred from the dispersion mea- 
sure is much larger than that of J1741— 2054. The 
search for radio counterparts to the LAT-detected pul- 
sars is still ongoing, the small number of radio detec- 
tions reported so far nevertheless indicates that a large 
fraction of LAT-detected pulsars are radio-quiet. 



3. What do we learn? 

To summarize the Fermi LAT observations of pul- 
sars during its first year of operation, a total of 55 
pulsars have so far been detected in gamma rays with 
high confidence. These include: 



quiet gamma-ray pulsars is an important discrimina- 
tor for theoretical models. The LAT discovered pul- 
sars are not necessarily radio-quiet, and thus need 
to be searched deeply for radio pulsations, helped by 
the pulsar positions and rotational parameters derived 
from the gamma-ray observations. So far, radio pul- 
sations have been detected for three of them: PSR 
J1741-2054 and J2032-H4127 at the Green Bank Tele- 
scope (GET) [13 and PSR J1907-h0602 at the Arecibo 
telescope Figure [5] shows the two-peaked gamma- 
ray light curve of PSR J2032+4127 as observed with 
the LAT, and the phase-aligned radio emission profile, 
providing the radio to gamma-ray phase separation. 



• 24 previously unknown normal pulsars, 

discovered in the LAT data. Among these 24 
gamma-ray pulsars, pulsations at radio wave- 
lengths have been detected for three. 

• 31 radio and X-ray loud pulsars, detected 
in gamma rays using rotational parameters mea- 
sured at radio telescopes. This number includes 
22 normal pulsars (among which are the six 
EGRET pulsars) and 9 millisecond pulsars. 

For both normal and millisecond pulsars, most 
gamma-ray light curves have two sharp peaks, while 



eConf C091122 



2009 Fermi Symposium, Washington, D.C., Nov. 2-5 



5 




Phase (|) 



Figure 6: Gamma-ray (upper panel) and radio (lower 
panel) light curves for the J0030-I-0451 millisecond 
pulsar, fitted with the Outer Gap (green) and Slot Gap 
(magenta) emission models. The observed emission is 
well reproduced by the two high-altitude emission 
models. Figure taken from [4l| ]. 



a few objects have one broad peak. The gamma- 
ray emission beam then seems to consist of a fan- 
like beam. Depending on the viewing geometry, this 
beam can either intersect the line of sight of the ob- 
server twice, resulting in two sharp peaks, or graze 
the line of sight, resulting in a broad structure. This 
emission geometry matches well the prediction of the 
Outer Gap (OG) and Slot Gap (SG) models, for which 
the gamma-ray emission is produced at high altitude 
in the magnetosphere. In these models, gamma-ray 
pulsars emit high-energy photons in fan-like beams, 
that cover large fractions of the sky. The evidence 
for high-altitude emission in OG or SG scenarios are 
supported by spectral measurements. The emitted 
spectra are consistent with exponentially cutoff power 
laws, with cutoff energies below 10 GeV in all cases. 
As mentioned above, this spectral shape reveals the 
absence of magnetic pair attenuation, expected in case 
of emission produced at low altitude, where the mag- 
netic field is most intense. Thus, both light curve and 
spectral properties favor outer magnetospheric emis- 
sion for gamma-ray pulsars. 

The observed similarities between normal and mil- 
lisecond pulsars in terms of gamma-ray emission in- 
dicate that the same mechanisms are operating for 
both pulsar classes. Indeed, the Outer Gap and Slot 
Gap emission scenarios fit the observed light curves 
of g amma-ray millisecond pulsars well in most cases 
[4lL |4^ . Figure [6] shows the observed radio and 
gamma-ray light curve for PSR J0030-I-0451, and fits 
obtained with the Outer Gap and Slot Gap emission 
models. These fits provide valuable constraints on the 
pulsar geometry (magnetic axis and line of sight in- 
clination angles relative to the rotation axis, a and 



C), generally in agreement with independent measure- 
ments, when available. Interestingly though, the PSR 
J1744-1134 and J2124-3358 MSPs are exclusively fit 
by the Pair Starved Polar Cap (PSPC) model, for 
which the gamma-ray emission is produced above the 
polar caps, at high altitude in the magnetosphere (see 
|4l|). High-altitude emission is thus favored for both 
millisecond and normal pulsars, although the location 
of the particle acceleration site may vary. 

After the observations of six pulsars in gamma rays 
with the EGRET telescope, there was evidence that 
the gamma-ray luminosity was proportional to 

\/e, where E is the rotational energy loss rate via 
electromagnetic braking [3^. The question of the 
conversion of spin-down energy loss rate into gamma- 
ray luminosity was therefore a key question to be ad- 
dressed with the LAT. Knowing the integral energy 
flux measured at the telescope, h, and the pulsar dis- 
tance d, one can calculate the gamma-ray luminosity 
Lj, given by: 

= iirfnhd^ (1) 

In this expression, /n is a fiux correction factor, 
which depends on the pulsar configuration (a and ^ 
angles) and the assumed beaming geometry [43| . This 
factor reflects the fact that the flux emitted towards 
the Earth may be different than that averaged over 
the sky, and is close to 1 for normal pulsars in outer 
magnetospheric emission models, with important vari- 
ations though (see Table 1 of [43|). This factor is also 
shown to be close to 1 for gamma-ray millisecond pul- 
sars [IJ ■ Figure [7] shows the gamma-ray luminosity 
for pulsars with known distances, assuming an /n of 
1 for all pulsars. The luminosity values suffer impor- 
tant uncertainties, because of badly known distances 
or emission geometries. Nevertheless, the gamma-ray 
luminosity trend seen with EGRET for high E pul- 
sars seems to stand. For the low E objects, includ- 
ing MSPs for which precise parallax distance measure- 
ments are generally available, the gamma-ray luminos- 
ity more seems to correlate with E rather than with 

\/e, indicating an inflection in the luminosity versus 
spin-down rate relation. Several objects in this low E 
regime are however found to be above the — E line, 
implying an absurd efficiency of conversion rj = L-y/ E 
greater than 100%. Their distances are probably over- 
estimated. More precise distance estimates, and more 
gamma-ray pulsars populating the 10'^'* to lO'^^ erg/s 
region of this plot might help see where the inflection 
occurs. 

As most known pulsars have spin-down energy loss 
rates below 10'^^ erg/s, it is reasonable to assume the 
Lj oc E relationship and thus use E/d'^ as a criterion 
for pulsar visibility; although it obviously lacks the a 
priori unknown beaming geometry information. Fig- 
ure [5] shows the "spin-down flux" E/d'^ as a function of 
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Figure 7: Gamma-ray luminosity above 100 MeV as a function of the spin-down energy loss rate E, for gamma-ray 
pulsars with known distances. Pulsars with two distance estimates have two markers connected with dashed error bars. 
Red triangles show the gamma-ray millisecond pulsars, blue squares show the LAT-discovered pulsars and green circles 
are normal radio-loud pulsars. Pulsars whose distance is based on the dispersion measure only are shown with open 
symbols. Additional details on this plot can be found in [l4| ]. 
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Figure 8: "Spin-down flux" E/d^ as a function of period for pulsars with known distances. Red triangles show the 
gamma-ray millisecond pulsars, blue squares show the LAT-discovered pulsars (accurate distances are unavailable for 
most of these pulsars, so that only a few of them are represented) and green circles are normal radio-loud objects. Gray 
dots are currently undetected pulsars. For these objects, DM distances were assumed, using the NE2QQ1 galactic 
electron density model [SB]. When possible, E values have been corrected for the Shklovskii effect [sH] which artificially 
increases the apparent spin-down rate P of pulsars with important transverse velocity. Globular cluster pulsars are not 
represented in this plot. 
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period for pulsars with known distances. As expected, 
the detected pulsars lie in the upper part of the plot. 
The lowest MSP on this graph is PSR J1614-2230, 
whose distance is clearly overestimated, leading to an 
gamma-ray efficiency of more than 100%, as can be 
seen in Figure [T] The kinematic Shklovskii contribu- 
tion [3^, which makes the apparent spin-down rate 
P larger than its intrinsic value, is therefore overesti- 
mated as well, and the actual E value must then be 
larger than we think it is. With both a smaller dis- 
tance and a larger E, this pulsar should be higher 
in the E/cP versus P plot. Bearing this in mind, 
we see that nearly all MSPs above a certain E/(P 
threshold have been detected with the Fermi LAT. 
Similarly, a good fraction of high E/cP normal pul- 
sars are detected. Also, most LAT-detected pulsars 
are not shown in this plot because of unknown dis- 
tances, and they are expected to populate the upper 
part of the plot. Most detected normal and millisec- 
ond pulsars are found to be above a common E/cP 
threshold, again indicating similarities in gamma-ray 
emission mechanisms. 

There are a number of possible reasons for the ex- 
istence of gray dots in the upper part of Figure [51 
some pulsars may be faint and therefore need addi- 
tional LAT data to be detectable (for instance, low 
galactic latitude pulsars, currently undetectable be- 
cause of the intense galactic diffuse background). Ac- 
cumulated photon counts will doubtless increase the 
number of detected gamma-ray pulsars. Some high 
E/(P pulsars may be bright gamma-ray pulsars, for 
which beams do not point towards the Earth. Addi- 
tionally, some pulsars may well not emit in gamma 
rays at all. Detailed geometrical studies of individ- 
ual objects might help understand the causes of non- 
detection. Also, it should be emphasized that DM- 
based distances are uncertain and some undetected 
high E/(P pulsars may be more distant, explaining 
their gamma-mutism. 

The gamma-ray detectability does not seem to cor- 
relate with the value of the magnetic field at the stel- 
lar surface. Instead, the pulsars detected with the 
LAT have high values of the magnetic field at the 

light cylinder, Blc « \J P / P^ oc \Je/P'^. Figure [9] 
shows the values of iJ^c as a function of the rotational 
period P for known pulsars. As we can see on this 
plot, both normal and millisecond pulsars detected in 
gamma rays are found in the upper part of their re- 
spective population. This again suggests that outer 
magnetospheric properties play an important role in 
gamma-ray emission for pulsars. For emission tak- 
ing place in the outer magnetosphere, it is expected 
that the electric field parallel to the magnetic field 
lines is proportional to -BlC: in turn implying that the 
cutoff energy of the emitted spectra depends on B^c 
(see e.g. [32|, |46|). Such a correlation might indeed 
be observed (see Figure 7 of [l3l)j although strongly 



relying on the highest measured values of the cutoff 
energy Ec^ which are rather uncertain. A reanalysis 
of these spectral parameters with increased datasets 
will reduce the uncertainties and make the relation- 
ship between cutoff energy and magnetic field at the 
light cylinder clearer. 

4. Search for off-pulse emission 

Some young and energetic pulsars like the Crab are 
known to convert a large portion of their spin-down 
energy loss rate into a wind of energetic particles, the 
so-called "Pulsar Wind Nebulae" (PWNe), producing 
synchrotron and inverse-Compton emission. A num- 
ber of these objects have been detected in the TeV 
domain, but the Crab Nebula was the only object 
ever detected in CeV gamma rays until Fermi was 
launched. The off-pulse emission of the Crab pulsar, 
that is, the continuous emission outside of the pul- 
sar gamma-ray peaks, has been analyzed using eight 
months of LAT data [IjI . The spectrum of the nebula 
is found to be well fit by the addition of a synchrotron 
and an inverse-Compton component, with no clear 
energy cutoff appearing for the synchrotron compo- 
nent. Figure [TUl shows the spectral energy distribution 
of the Crab Nebula from soft to high-energy gamma 
rays. The spectrum measured with the LAT connects 
with lower and higher energy results within error bars, 
allowing for cross-calibration with Cherenkov experi- 
ments. 

Besides the observation of the Crab Nebula with 
the LAT, almost half of the pulsars detected during 
the first year are found to be spatially associated with 
known X-ray or TeV PWNe and Supernova Remnants 
(SNRs) [13, dl- Some gamma-ray pulsars are asso- 
ciated with TeV sources not having counterparts at 
other wavelengths, which may point to yet unknown 
PWNe or SNRs. For instance, the association of the 
Fermi-detected pulsar PSR J1907-h0602 with the TeV 
source MGRO J1908-I-06 suggests that the latter is 
the pulsar wind nebula of the former [l^ . A system- 
atic search for continuous (non-pulsed) emission in the 
off-pulse of gamma-ray pulsars is underway, and may 
foster detections of new PWNe or SNRs. Significant 
gamma-ray emission is already seen from objects other 
than the Crab nebula: the LAT detected GeV emis- 
sion from the Vela-X pulsar wind nebula [l7[ , provid- 
ing constraints on the injected populations of charged 
particles. Another example is the detection of GeV 
gamma-ray emission from the pulsar wind nebula in 
MSH 15-5Zj)owered by the COMPTEL pulsar PSR 
B1509— 58 [23]. FigurefTTlshows a map of the observed 
emission above 10 GeV in one year of data. An excess 
is clearly observed around the PSR B1509— 58 pul- 
sar, coincident with H.E.S.S. contours of the pulsar 
wind nebula. Other such detections of nebular emis- 
sion may be reported in the next months, and Fermi 
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Figure 9: Magnetic field at the light cylinder Blc as a function of the rotational period P for known pulsars. The 
symbol code is the same as in Figure [5] When possible, the kinematic Shklovskii effect [s^ has been taken into account 
in the calculation of Blc values. Globular cluster pulsars are not represented in this plot. 




Figure 10: Spectral energy distribution of the Crab Nebula, along with results of other experiments from soft to very 
high-energy gamma rays. See [l3| for details on the analysis, on the fits of the inverse Compton and synchrotron 
emission components and for references. 



LAT data along with multi-wavelength observations 
might help understand the physics of these objects 
associated with energetic pulsars. 



5. Conclusion and prospects 

Pulsed gamma-ray emission has been observed for 
55 pulsars over the first year, including previously un- 



known pulsars and millisecond pulsars. These pulsars 
are generally close and energetic (high E/cP values), 
and have high magnetic fields at the light cylinder, 
Bi^c- Their gamma-ray luminosity increases with 
the spin-down energy loss rate E. Low E pulsars seem 
to follow oc E, whereas high E ones rather fol- 
low oc \/e. Refined spectral measurements and 
more pulsar detections will help see where the inflec- 
tion occurs. Light curve shapes and spectra are com- 
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Figure 11: Gamma-ray counts map around the MSH 15—52 pulsar wind nebula above 10 GeV, with H.E.S.S. contours 
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parable for normal and millisecond pulsars, indicating 
that the same mechanisms are operating. Gamma-ray 
pulsars seem to emit fan-like beams, separated from 
the radio emission (if any) , and produced in the outer 
magnetosphere. Outer Gap and Slot Gap models are 
therefore favored over the Polar Cap model, although 
more data are required to discriminate between the 
two outer magnetospheric emission models. In addi- 
tion, a few millisecond pulsars seem not to follow the 
Outer Gap and Slot Gap models, but are rather fit 
with a third family of high-altitude emission models, 
the Pair Starved Polar Cap model. Much has already 
been learned about the way gamma-ray pulsars emit; 
but unexpected features may appear in future pulsar 
observations, requiring refinement of current models. 

The Fermi LAT was initially expected to operate 
for five to ten years. If this goal is attained, then the 
LAT will eventually be two to three times more sensi- 
tive than it was during the first year. As can be seen in 
FigurelHl the LAT has so far detected the highest E/(P 
objects, and a lot of distant and/or low E pulsars 
may be awaiting detection, provided their gamma-ray 
beams are pointing towards the Earth. In addition to 
the simple accumulation of gamma-ray data, another 
way to increase the LAT sensitivity is to improve pul- 
sation search algorithms. Such analysis refinements 
are being developed, increasing the sensitivity of the 
LAT at low energies or enhancing the pulsation search 
by giving more weight tophotons that are likely asso- 
ciated with pulsars [U [30] ■ 

More than half known MSPs are in globular clusters 
and are not represented in Figure[Sl None of these pul- 



sars have yet been detected in gamma rays with the 
LAT, most probably because of their large distances. 
However, the LAT has detected continuous emission 
from the close globular cluster 47 Tucanae 3 (see Fig- 
urc [T2|) . The observed spectra can be explained by the 
cumulative emission of seven to sixty- two MSPs. The 
LAT might eventually be able to see individual pul- 
sars in globular clusters. Also, the blind searches of 
gamma-ray pulsars in LAT data have not been search- 
ing for unknown MSPs thus far. Some of these cur- 
rently unknown MSPs may be radio-faint, and rapidly 
rotating radio-faint pulsars have never been detected, 
although they have all reasons to exist. 

Finally, there is much to expect from the follow- 
up observations of LAT unidentified sources at other 
wavelengths, such as those of the Fermi LAT Bright 
Sources List First searches for pulsations at 

the position of LAT unidentified sources have already 
been extremely successful in finding new pulsars: for 
instance, several MSPs have been found at radio wave- 
lengths with the Green Bank Telescope [s^. Similar 
searches are being conducted at Effelsberg, Nangay, 
Jodrell Bank, Arecibo and Parkes radio telescopes. 
The discovered pulsars can then be monitored at radio 
wavelengths, and searched for gamma-ray pulsations 
in LAT data afterwards, using accurate rotation pa- 
rameters. Radio telescopes have been leading pulsar 
science over the last 40 years, and the pulsars detected 
at optical. X-ray or gamma-ray wavelengths are still 
a minority compared to the '^2000 observed in ra- 
dio. For the first time though, LAT is providing radio 
telescopes with a nearly unbiased way to search for 
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Figure 12: 1.5° x 1.5° region around the 47 Tucanae globular cluster, as seen with the LAT in eight months of data. 
Black contours show the stellar density in the cluster, while the white circle indicates the 95% confidence region of the 
gamma- ray source. Additional details on this image can be found in Q]. 



energetic normal and millisecond pulsars, by guiding 
them to unidentified gamma-ray sources in all direc- 
tions of the sky. In addition to helping understand 
the physics of gamma-ray emission in pulsar magneto- 
spheres, the LAT is also helping discover new pulsars 
at other wavelengths, some of them being of poten- 
tial interest for gravitational wave detection projects, 
such as LEApi or NANOGrav^ or any project in- 
volving pulsars. We can therefore expect the Fermi 
LAT to provide an important contribution to pulsar 
astronomy during the next few years. 
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